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ABSTRACT 

Beam radii for cone-dominant pulsars foUow a power-law relation with frequency, d = {v / V())'' + -do, which 
has not yet well explained in previous works. We study this frequency dependence of beam radius (FDB) for 
cone-dominant pulsars by using the curvature radiation mechanism. Considering various density and energy 
distributions of particles in the pulsar open field line region, we numerically simulate the emission intensity 
distribution across emission height and rotation phase, and get integrated profiles at different frequencies and 
obtain the FDB curves. For the density model of a conal-like distribution, the simulated profiles always shrink 
to one component at high frequencies. In the density model with two separated density patches, the profiles 
always have two distinct components, and the power-law indices k are found to be in the range from -0.1 to -2.5, 
consistent with observational results. Energy distributions of streaming particles have significant influence on the 
frequency-dependence behavior Radial energy decay of particles are necessary to get proper in models. We 
conclude that by using the curvature radiation mechanism, the observed frequency dependence of beam radius for 
the cone-dominant pulsars can only be explained by the emission model of particles in two density patches with a 
Gaussian energy distribution and a radial energy loss. 

Subject headings: pulsars: general - star; magnetic fields - relativistic particles - curvature radiation 



1. INTRODUCTION 

A pulsar profile is explained as being the line-of-sight cuts 
through the emission beam when a pulsar rotates. For cone- 
dominant pulsars, profiles widen at low frequencies. The varia- 
tion of pulsar profile width or component separation, W , against 
the obser vation frequenc y, v, can be described by a power-law 
function dThorsettll 19911) 



Table 1 

Geometry parameters and fitted k and -do for 7 

CONE-DOMINANT PU LSARS WITH CLEAR FREQ UENCY 
DEPENDENCE dMlTRA & RANKINl2002b . 



W = av'' + WQ. 



(1) 



Here, Wo is a width constant. With a dipole field geometry 
for pulsar emission region, the emission of lower frequencies 
is believed to be generated at larger heights for the wider open 
angles of the emission cone. This is well known as the radius- 
to-frequency mapping (iRuderman & Sutherlandl[r975l: ICordesI 
[1978; Phillips 1992). 

Observed profile width can be expressed by, for example, the 
peak separation between the outermost components, Wpp (i.e. 
the peak-peak separation), the pulse width at the 10% of the 
highest intensity peak, Wio, and the pulse width at the 50% 
of the peak, W50. In order to understand the frequency de- 
pendence of pulsar p rofiles in a much wider frequency range, 
IXilouris et al.l (1 19961) observed a number of nearby bright pul- 
sars up to 32 GHz and measured their W50. Combining with pre- 
vious measurements at lower frequencies, they fitted the pulse 
widths by W50 = au'^+Wo, and got -0.29 >k> -0.94 for a sam- 
ple of pulsars. 

The radius of the pulsar emission beam, d, is related to the 
profile width, W, by (Gil 1981), 
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1? = 2arcsin[sin^ sinQ!sin(a + /3) + sin'^ 



(2) 



Here, a is the inclination angle of the magnetic axis related to 
the rotation axis, /3 is the impact angle of the line of sight to 
the magnetic axis. For many pulsars with good polarization 
measurements, the values of a and /3 have been determined 
jLvne & Manchester [19881; iRankinI IT990t lEverett & WeisbergI 
120011) . For a given pulsar, the beam radius, and the pro- 
file width, W , are quasi-linearly related for various frequencies. 



Therefore, the variation of pulse-width with frequency should 
be physicall y related to the frequen cy dependence of pulsar 
beam radius. iMitra & RankinI (l2002h collected observed pulse 
widths for 7 cone-dominant pulsars and calculated their beam 
radii at various frequencies. They found that the radius of the 
outer-cone beam of 7 pulsars follows 

d = (vlvof + do, (3) 

here t?o is the beam radius at the infinite-frequency, and the vq is 
the characteristic frequency. The values of the power-law index 
k are in the range of -0.3 to -1.2. The values of (3 are smaller 
than those of do, as shown in Table [T] We noticed that for a 
given pulsar the k value obtained from pulse widths is almost 
the same as that fitted from the beam radii. iMitra & RankinI 
(2002) noticed that the beam radii of inner-cone do not show 
the frequency dependence. 

The frequency dependence of beam radii or profile widths 
can be easily explained by the open field-lines of the emission 
region in a dipole field geometry of neutron stars, as lon g as 
the radius-frequency-mapping holds (Ruderman & Suthe rland 
IT975I) . Different emission mechanisms can lead to different 
power-law indices k. The electron bremsstrahlung model (IVitarmo & Jauh 
^73) predicts k ^ -0.45. The vacuum inner gap model (Ruderman & SutI 
T975) can give A: ~ -1/3. The curvature plasma model dBeskin et al.l 
il988i) can have two extremes for the ordinary emission mode, 
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Fig. 1. — Geometry and parameters defined for an emitting beam. Vindi- 
cates the rotation axis in z direction, /i represents tlie magnetic axis, which is 
inclined by an angle of a with respect to Q. The LOF is tangential to the light 
cylinder, which has a polar angle 9 from the magnetic axis /i. The radiation 
beam has a radius of i? at a height of r. 



either k ^ -0. 14 or A: ^ -0.29. T he emission from the cyclotron 
instabiUty dMachabeli & Usovl ll989) gives k ^ -0.17. These 
theoretical values can not cover the wide range of the observed 
k. 

In this paper, we try to explain the frequency dependence of 
the radio beam observed for cone-dominant pulsars using the 
curvature radiation mechanism. If the edge of pulsar beam is 
generated by particles flowing along the last open field-lines 
(LOF), we can calculate the radio beams, and investigate their 
dependence on pulsar parameters, period P, inclination angle 
a, impact angle /3 and the Lorentz factor of particles 7. Using 
some generalized energy and density distributions of particles 
in the magnetosphere, we numerically calculate the radio emis- 
sion beams and fit their frequency dependence. We also inves- 
tigate the influence on the pulsar beam by the radial decrease of 
particle energy and the particle energy distribution. 

2. CURVATURE RADIATION BEAM FOR PARTICLES WITH A 
GIVEN 7 IN A DIPOLE FIELD 

In general, pulsar radio emission is assumed to be generated 
by curvature radiation of secondary particles streaming along 
the last open field-lines. In the radio emission region, mag- 
netic fields can be described as a static dipole (e.g. lGangadharal 
I2OO4I) . because the multipolar field components of a neutron 
star vanish there and the swee ping effect due to rotation is also 
negligible according to lDvks & HardingI (|2004|) . Therefore, in 
this paper, we will use the dipole field to study the frequency 
dependence of beam radius (FDB). 

The size and geometry of a dipole magnetic field is deter- 
mined by pulsar period P, inclination angle a and impact angle 
p. In the polar coordinate system with the polar axis along the 
magnetic axis direction (Fig. [U, a dipole field-line can be de- 
scribed by r = r^, sin^ 9, here 6 is the polar angle from the mag- 
netic axis, r is the distance from the dipole origin, re is the field- 
line constant, which is the distance from the origin to the point 
of the field-line intersection with the magnetic equatorial plane 
of = 90°. For an inclined dipole, the LOF are contained in 
the light cylinder (see Fig.lTJ. The radius of the light-cylinder, 
7?ic = cP/2ti, gives the limit of the field line constants rg for the 



Fig. 2. — Frequency dependence of the beam radius (FDB curve) calculated 
from the curvature radiation from particles of a single 7 for an inclined dipole 
of a = 60°. Beam radii in the magnetic azimuth of <j> = 0° and (p = 90° are 
frequency dependent, though the beam size of = 0° is smaller than that of 
(p = 90° at any frequency. This behavior is shown in the top panel for 300 
MHz, 1.0 GHz and 2.0 GHz emission. The solid lines are calculated using two 
terms in Eq. Q, and the dashed lines are only for its first term. Here period 
P = 1 s and 7 = 400 are taken for calculations. 



LOF, which is different for pulsars with different periods. The 
angular diameter of the polar cap defined by the feet of LOF 
on the neutron star surface is related to the pulsar period, P, by 
20pc = 1.6°P~'/^. The opening angle of the beam from the tan- 
gents of the LOF near the surface is about 1 .2°f which de- 
fines the minimum geometrical beam angle. Observations show 
t?o > 1.2°p-i/2. The emission beam determined by the LOF is 
not circular but compressed in the meridional direction in the 
plane of rotation and magnetic axes (Biggs 1990). We define 
the magnetic azimuthal angle, cf), starting from the connection 
between the magnetic axis and the rotation axis (to the top be- 
ing the north) as being (p = 0° (see Fig. |2]i. The beam radius in 
the direction of = 0° is smaller than that in the direction of 
= 90° (to the east). 

In the curvature radiation mechanism, the emission frequency 
is not only related to the field geometry but also to the Lorentz 
factor 7 of particles. The simplest case we consider here is 
that particles have the same Lorentz factor 7, and that the radio 
beam is defined by the tangents of the LOF. When a relativis- 
tic particle streams along a field-line, it can produce curvature 
radiation with a characteristic frequency of 

- ^^^^ (4) 
47rp 

Here, 7 is the Lorentz factor of particles in the range 10^- 10"*, 
p is the curvature radius of the particle trajec tory. In any field- 
line, the curvature radius can be expressed by (IGan gadharatZ 0041) 
sin6i(5 + 3cos26')3/2 

3\/2(3 + 2cos26') 
Note that 9 varies with r. The angle between the tangent of a 
field-line of 9 and the magnetic axis is 



, l + 3cos2e' , 

V = arrrn<i( ) 



(6) 



VlO + 6cos26i 

Combining Eq. ^ and Eq. one can find the relation be- 
tween z9 and p for any field-line of r^. Because the radiation 
frequency i/ is related to p by Eq. (|4|l, we get 

3 3 
^= 15.4°(2^) + 0.43°(^)l (7) 
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Fig. 3. — Beam radii at a series of frequencies calculated from the curvature 
radiation of particles with a single 7. Default model parameters are P = 0.5 s, 
a = 60° , /3 = 0° and 7 = 400. Various curves have been calculated for different 
parameters, a, fi, P and 7 in four panels. Note that the two curves for a = 60° 
and a = 30° in the upper left panel are completely overlapped. 



One can use Eq. (15) from 'Gangad haral (|2004|) to get the field- 
line constants for each LOF. Obviously, as indicated in Eq. (|7]l, 
the beam radius is related to the radiation frequency. At a higher 
frequency of ^ 'y^c/re, the second term can be neglected, 
and the beam radius ■& is related to the frequency ly by roughly 
a power-law with the index of k = -I. At a low frequency of 
v < 7^c / r^, the frequency dependence becomes slightly steeper 
due to the contribution from the second term. Note that the 
beam radii in all magnetic azimuthal directions are frequency 
dependent. Clearly, the first term of Eq. (|7} is a good approxi- 
mation of the frequency dependence of pulsar beam, which has 
a power-law index k = -l. 

We calculated the beam radii at different frequencies for var- 
ious model parameters, as shown in Fig.|3] All curves show the 
frequency dependence of radiation beam with a power-law in- 
dex of approximately -1, which is diffe rent from observationa l 
values in the range of -0.3 >k>-l.2 dMitra & Ranldnll2002l) . 
We found that there is no influence on the frequency depen- 
dence of beam size by the magnetic inclination a, because a 
only leads to the compression of beam in the meridional direc- 
tion, and has almost no influence on the radius of 4> = 90°. Dif- 
ferent impacts of line of sight with different /3 values lead to dif- 
ferent cuts of the beam from /3 = 0° to /3 = The lower limit of 
1} is determined by the angular size of the polar cap {3/2dpc) or 
the P value. Pulsars with a small period P have a small R\c and 
hence a small and a small curvature radius [Eq. (|5)], which 
corresponds to a larger radiation frequency [Eq. (3]], as shown 
in Fig. |3] Note also that in the pulsar emission region curva- 
ture radiation of particles with a single 7 cannot produce radio 
emission over a wide frequency range from hundreds MHz to 
ten GHz. Particles with larger 7 produce curvature radiation at 
much higher frequency, which shifts the FDB curves to higher 
frequency ranges. 

The calculations shown in Fig|2] and |3] were made with as- 
sum ptions that pulsar beams are bo unded by the LOF. How- 
ever, |^deman]&^utherlan2([l975l) suggested that beam edge 
should be bounded by the critical field-lines, which are orthogo- 
nal to (instead of tangent to) the light cylinder at the intersection 
points. The critical field-lines are located between the magnetic 
axis and the LOF. Here, the parameter ry is used to describe the 
location of field-lines, with ry = for the magnetic axis, 77 = 1 
for the LOF. To generate the curvature emission beam of the 



same open angle, the emission height from the critical field 
lines (0 < ?7 < 0.74, depending on a) is larger than that from 
the LOF as shown in Fig. |3] and the curvature radius is also 
larger, so that the emission has a smaller frequency of, 

373c 



z/ = X rj 



A'Kp 



(8) 



with i^iof the corresponding emission frequency for the LOF. We 
carried out a set of calculations, and found that the curvature ra- 
diation for the critical field lines has almost the same frequency 
dependence of emission beam as that for the LOF, and k and 
do values are consistent within 5%, though i/q is smaller due to 
larger curvature radii of field lines. 

3. CURVATURE RADIATION BEAM FROM PARTICLES WITH 
VARIOUS ENERGY AND DENSITY DISTRIBUTIONS 

Particles in pulsar magnetosphere should have an energy dis- 
tribution, which radiate in a range of frequencies at a range of 
heights for a given LOF. Furthermore, particles flow out along a 
set of open field-lines, rather than just the LOF. The pulsar radio 
emission from a given height and a given rotation phase is con- 
tributed from particles not only in the field lines which are tan- 
gential towards the observer, but also in the nearby field lines in 
the bunch within the 1 /7 emission cone. The emission is coher- 
ent ra diation from a bunch of particles (Buschauer & Benfordj 
11971) . 

According to simulation results of iMedin & Lail (1201 Ol) . we 
assume in this section that secondary particles for curvature 
radiation at radio bands follow a Gaussian energy distribution 
with a peak at 7,„ of several hundreds: 

ne(7)^exp[-^^^— ^]. (9) 

Here, the standard deviation cr^ is of several tens. Considering 
the continuity of the particles flowing along the field-line tube, 
we got the number density of particles at r as 

n,(r) = n,o(r/R.)-\ (10) 

Here, n^o represents the number density at the bottom of a mag- 
netic field tube near the surface of a neutron star. 

The power of curvature radiation at a frequency ly from one 
particle is given by. 
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-i-r 

p 



(11) 



Nb particles in a field bunch in the region with a dimension size 
of less than half emission wavelength produce the total emis- 
sion power of approximately NI P^. 

Because the curvature radius varies everywhere in the dipole 
field, according to Eq. (|4]i and (fTTT i. the observed emission at 
frequency u and rotation phase (p should come from the tan- 
gents of a set of open field-lines with the same magnetic az- 
imuth (f), where the curvature radius p of a field-line and par- 
ticle energy 7 are nicely combined to produce emission. The 
observed total power therefore should be the sum of them. 



n^(r,j)NiP^ir,(l),-f)dr. 



(12) 



Note that the magnetic azimuth is a mo notonic function of 
the rotation phase ip, as given by Eq. (11) in lGangadharal (|2004|) 
which is equivalent to the well-known 5-curve for polarization 
angle. 
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Fig. 4. — Illustration for the uniform density distribution of particles in 
pulsar magnetosphere. Here, = n^Q(r/Ri,y^. 
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(a) without A&R effects 



(b) with A&R effects 



Fig. 5. — (a): Upper panel: Emission intensities at a given frequency (here 
V = 200 MHz) for various rotation phases ip and emission heights, radiated by 
a group of particles which have a Gaussian energy distribution and an uniform 
number distribution in magnetic polar angles as shown in Fig.O Lower panel: 
The accumulated emission profile (also al v = 200 MHz). The model param- 
eters used are /" = Is, a = 45° and /3 = 3°, and 7,,, = 350 and cr^ = 30. (h): 
Same as left panels (a), but the abberation and retardation effects have been 
considered. 



3.1. Emission beam from particles uniformly distributed in 
magnetic polar angles 

As the first try, we assume that particles follow a Gaussian 
energy distribution and are uniformly distributed in the open 
field-lines for all magnetic polar angles, as shown in Fig. |4] 

We can calculate the emission from each open field-line at 
every tangent for each frequency. At a given emission height, 
the central open field-lines have larger curvature radii than the 
outer open field-lines. Particles with a given Lorentz factor in 
the central field lines would produce a lower frequency emis- 
sion [see Eq. (|4]l] than they do in the outer field lines. For 
a group of particles with a Gaussian energy distribution, only 
these particles with larger Lorentz factors in the central field 
lines can produce the emission at the same frequency as parti- 
cles of lower Lorentz factors in the outer field lines, but their 
emission intensity is reduced by a factor of 1 /7^. Note that 
magnetic fields have smaller curvature radii at lower height. All 
these factors produce the bifurcation feature shown in FigJS] 

We limited our calculations for the region of r^ax < 0.1/?ic, 
because the emission intensity from higher regions is negligible 
(see Fig. Is) due to the relatively small particle density and large 
curvature radius. For P = \ s, /?ic ~ 47771 km c:; AlllR^,, if 
we ta ke iJ* =10 km. We consider the phase shift (iGangadharal 
l2005h of the emission from any height of every field-line, and 
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magnetosphere. The density follows Eq 



integrate their emission for profiles at every observation fre- 
quency (lower panels in Fig. |5]l. We found that aberration and 
retardation effects can change the projection position of emis- 
sion in the sky plane, which cause the profile to be enhanced 
by about 20% at the phase of = -10° and weakened by about 
22% at the phase of (/j = 10° for the model calculations in Fig.|5] 
We conclude that curvature radiation from particles uniformly 
distributed in magnetic polar angles cannot produce double conal 
profiles. 

3.2. Emission beam from particles in the conal density 
distribution 

Small curvature radii of the magnetic field lines near the 
edge of the pulsar polar cap make curvature radiation more ef- 
ficiently being generated, which then causes the sparkings more 

probably located around the edge of polar cap (Ruderman & Sutherland! 
'1975VOiao et al."2001). It is therefore possible that accelerated 
particles are more likely distributed in a conal area, instead of 
an uniform distribution. Noticed that charged particles flow out 
along a fixed line tube, we define the density distribution in a 
conal area on the polar cap as 

neo(??) = neoexp[ ] , (13) 

here r]„ denotes the peak position of the cone, and the standard 
derivation ct^ describes the cone width. At any radius r and any 
magnetic polar angle rj the density is 

n,(r,r]) = n^(r])(r/R,r\ (14) 

which does not vary with magnetic azimuthal angle (j). The 
illustration for the conal density distribution is given in Fig.|6] 
Similar to the uniform model in Sect. 3.1, we calculated the 
curvature radiation of particles in such a density cone from each 
field line at every tangent at each frequency, to composite the 
emission beam between 100 MHz to 30GHz. Fig.|2]shows the 
distribution of emission intensities for various heights and ro- 
tation phases as well as integrated profiles at two example fre- 
quencies, 200 MHz and 800MHz. There are two components 
in the integrated profile at 200 MHz, though the emission for 
central rotation phases coming from the region of several tens 
of Ri, is relatively strong compared to those from higher re- 
gions. The aberration and retardation effects, which shift high- 
altitude emission to early rotation phases, can cause the lead- 
ing profile component to be wider and weaker than the t rail- 
ing one at 200 MHz, as demonstrated bv lDvks etall (|2010'). At 
800 MHz emission mostly comes from a lower region of about 
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(a) 200 MHz 



(b) 800 MHz 



Fig. 7. — Same as Fig. |5]for emission intensity distributions and tlie accumu- 
lated profiles (a): at 200 MHz and (b): at 800 MHz for a conal density model. 
The particle energy and density distribution parameters are 7,,, = 350, cr^ = 30, 
ri,„ = 0.75 and cr^ = 0.08, with P = Is, a = 45° and /3 = 5°. 




Fig. 8. — Integrated profiles at a series of frequencies calculated from curva- 
ture radiation of particles with a conal density distribution. The model param- 
eters are the same as those in Fig. [7] The profiles are aUgned to the absolute 
zero rotation phase when the sight line is in the meiidional plane defined by 
the rotation and magnetic axes. 



r = 10 - 2QR^, rather than a few hundred of for 200 MHz 
emission, where the sight line cuts only the beam edge and the 
emission from higher regions is too weak, so that only one com- 
ponent appears in the integrated profile. The profile evolves 
from two components at low frequencies to one component at 
frequencies of above 400 MHz (see Fig.[8]l. 

We measured the profile width at 10% of the peak intensity, 
and calculate the corresponding beam radius -d using Eq. (|2]i. 
The FDB curves are plotted in Fig. |9]for various initial con- 
ditions. These curves are fitted using Eq. (|3]l, and results are 
listed in Table|2] The profile width is close to zero when /3 = i?. 
We have to fix i?o to be the same as the impact angle /3 dur- 
ing the fitting, because the emission at very high frequencies 
comes from the lowest emission point on the LOF. As seen in 
Fig.|9]and Table. |2] the FDB curves are closely related to emis- 
sion geometry (a, particle energy (7, a-y) and conal density 
distributions (rj,,,). 

We noticed that the FDB curves in all panels of Fig. |9] are 
converged to a similarly flattened FDB curve at high frequen- 
cies, except for various (3 in panel b. As leaned from Fig.|7] the 
emission at high frequency comes from a very small height, so 
that the emission beam is rather narrow due to the geometrical 
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Fig. 9. — Various FDB curves calculated from particles with a Gaussian 
energy distribution and a conal density distribution above the polar cap with 



different model parameters of a, /3, 7^ 
ters ai'e the same as those in Fig.|2] 



(Tj , ri„ and . The default parame- 



TABLE2 

Fitted parameters for various FDB curves in Fig. [9] 

FROM VARIOUS CONAL DENSITY MODELS WITH DIFFERENT 
PARAMETERS OF a, /3, 7„,, a^, rj,,, AND a,, 
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effect. The density distribution {rjm, cr-q) and energy distribution 
(7, fj-y) of particles as well as the magnetic inclination (a) do 
not have any obvious effect on the beam at high frequencies, as 
shown by the similarly converged flat FDB curves in all panels 
in Fig. |9] As expected, a sight line with a small impact angle 
P can look into the much deeper magnetosphere and detect a 
much smaller beam at high frequencies (panel b in Fig.|9]l. 

The FDB curves in Fig.|9]are influenced by a, 7, and 77,,, 
more obviously at the lower frequency end. For a larger mag- 
netic inclination angle a, the polar cap as well as pulsar beam 
will be more compressed in the meridian dimension (Fig. |2). 
When other model parameters are fixed, for a dipole field ge- 
ometry with a larger a, the fixed 77,,, means that the bunch of 
field lines closer to the magnetic axis are used to calculate the 
emission beam, which have large curvature radii and emission 
frequencies are systematically smaller The characteristic fre- 
quency I'o of the FDB curves is therefore also smaller (see Ta- 
ble |2]l. Particles with a larger 7,„ produce emission of a given 
frequency at higher altitudes with larger curvature radii, which 
produce a wider beams (panel c in Fig. |9]l with a large char- 
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Fig. 10. — Distributions of k and ud from 6375 combinations of various 
model parameter values of the conal density model. 
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Fig. 1 1 . — Profile evolution with frequency for two extreme k values: Panel 
(a) for an extreme steep frequency dependence of k = —2.64 with model pa- 
rameters of a = 45°, /3 = 9°, 7,,, = 300, cr-y = 10, r),„ = 0.9 and ar, = 0.05; and 
Panel (b) for an extreme flat frequency dependence of k = —0.32 with model 
parameters of a = 65°, /3 = 1°, 7m = 200, a-y = 50, r]„ = 0.3 and cr,, = 0.25. 



acteristic frequency vq (Table |2]l. If the energy distribution of 
particles is rather wide, e.g. for a model with a larger a^, the 
emission region at a given frequency is very extended along 
field lines. However the dominant emission always comes from 
a lower emission height, and the integrated profile and hence the 
beam is much less sensitive to frequencies (panel d in Fig. |9]l. 
The particle distribution (77,,,) is directly related to the field lines 
where the emission is produced. It is easy to understand that 
particle distributed with a small 77,,, tends to emit with a smaller 
beam at any frequencies (panel e in Fig.|9]l. We noticed that the 
particle distribution width a.^ does not have obvious influence 
on the FDB curves ( panel/in Fig.|9]l. 

To further investigate the frequency dependence of pulsar 
beam radius on model parameters in a wide range of values, 
we have calculated the emission intensity distribution and inte- 
grated profiles (as shown in Fig. |7]and Fig. |8]l for 6375 combi- 
nations of different parameter values: a = (5°, 25°, 45°, 65°, 
85°), /3 = (1°, 3°, 5°, 7°, 9°), 7,„ = (200, 250, 300, 250, 400), 
cr^ = (10, 30, 50), Ju = (0.3, 0.5, 0.7, 0.9) and = (0.05, 0.1, 
0.15, 0.2, 0.25). Pulsar period is fixed to be P = Is. The fre- 
quency dependence of beam radii are fitted with Eq. ^ to get 
k and fQ. We got the distributions of k and i^o from these com- 
binations as shown in Fig. [TO] We found that values of k are 
in the range of -0.2 to -2.7, consistent with observed values. 
However, profile components always merge to one component 
above ^ 300MHz, which is very different f rom obser vational 
facts of cone-dominant pulsars ( Mitra & jRankinl l2002b . We 
show in Fig. [TT] the profile evolutions for two extremes of k 
values: k = -2.64 and k = -0.32. For the case of A: = -0.32, the 




Fig. 12. — Illustratio n fo r the particle distribution in two density patches. 
The density follows Eq. )15t . 
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Fig. 13. — Same as Fig. |5]for emission intensity distributions and the ac- 
cumulated profiles (a): at 400 MHz and (b): at 1.6 GHz for a model of two 
density patches. The particle energy and density distribution parameters are 
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= 350, 
Is, a = 



= 45' 



= 30, 
and ; 



: 0.75, 



: 0.08, 4,,,, = 45° and cr^ = 12.5°, with 



two profile components merge even at a very low frequency be- 
low lOOMHz. The merging happens when 7,,, < 250 or > 30 
or rjii, < 0.5 or (T,, > 0.15. We noticed that when 7,„ > 350, for 
models of curvature radiation from particles in the conal density 
distribution, the emission below about lOOMHz comes from the 
region of several thousands of 7?,^, which leads profiles as wide 
as about 200°. 

In summary, though curvature radiation models of particles 
in a conal density distribution can produce a proper range of k, 
the merging of profile components above few hundreds MHz is 
the main reason to reject the models for cone-dominated pro- 
files. 

3.3. Emission beam from particles of two density patches 

In the inner gap model (iRuderman & SutherlandllT975l) . the 
sparkings do not continually happen above the polar cap and 
are not concentrating in a ring area to produce the conal density 
distribution discussed above. The most probably case is that the 
sparkings occur separately above the polar ca p, which prod uce 
the subbeams o f the observed conal beam (iRankinI 1983h or 
the patchy beam (Lv ne & Manchested 19881: [Han & Manchested 
2001 ). In the following we work on the models of curvature ra- 
diation of particles in two density patches which could be rea- 
sonably produced by two sparkings in the inner gap of a rotating 
neutron star 

The density distribution of two patches can be defined on the 
pulsar polar cap region and extended to high magnetosphere. 
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Table 3 

Parameters fitted to the FDB curves in Fig.[T5]for models of two density patches with various a, j3, 7,,,, a^, rj„,, cr^, (j>, 
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Fig. 14. — Same as Fig.[8]except for the mo del o f two density patches. The 
model parameters are the same as those in Fig. 1131 



For simplicity, we model two density patches (see Fig.[T2]l which 
are symmetrically located about the meridional plane and peaked 
at (rim, (pm) and (77,,,, -0,,,), respectively, with Gaussian widths 
of (T,, and in the magnetic polar and azimuthal directions, 
i.e. 



"eC'-,'?, (/>) = neoir/R^) ^exp[- 



(ri-r]m) 
2^2 



]exp[- 



2al 



(15) 

The distribution of the emission intensity across height and 
the rotation phase and the evolution of integrated profiles with 
frequency for the curvature radiation model from particles of 
two density patches are shown in Fig. [13] and [14] respectively. 
Our model calculations show that: 1) the leading component 
is always slightly wider and weaker than the trailing one; 2) 
emission at a lower frequency comes from relatively higher and 
wider region (10- 50/?* for 400 MHz compared to 8 - 20/?* for 
1.6 GHz); 3) the low-frequency profiles become wider and are 
shifted towards an early rotation phase; 4) the two components 
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Fig. 15. — Same as Fig.|2]but for the mod el fo r two density patches. The 
model parameters are the same as those in Fig. 1131 



always remain resolved even at very high frequencies, because 
two density patches keep separated even down to the star sur- 
face. 

From simulated profiles at a set of frequencies, we calculated 
the beam radius and checked the frequency dependence of pul- 
sar beam. We also checked if FDB curves are related to emis- 
sion geometry (a, /3), particle energy distribution (7, a-,) and 
patch geometry (77,,,, cr,,, 4>„ and a^). The curves are shown in 
Fig.[T5]and fitted parameters are listed in Table|3] The influence 
of these parameters on the FDB curves are more or less simi- 
lar to those found from the conal density models. We noticed 
that (pm has the similar effect on the FDB curves as r]m, both of 
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which determine the separation of the two patches. Particles in 
two patches with a larger can produce the emission beams 
with a steeper FDB curve. 







Fig. 16. — Distributions of k and from 29070 combinations of various 
model parameter values of the two density patch model. 
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Fig. 17. — Profile evolutions with frequency in the two density patch model 
for two extreme values of k. Panel (a) for k = -2.02 and Panel (b) for k = —0.09. 
Model parameters for panel (a) are a = 45°, /3 = 5°, -fm = 400, (7^ = 10, r],,, = 
0.9, cTr, = 0.1, 0,„ = 65° and cr^ = 15°, and for panel (b) are o = 45°, ^ = 1°, 
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= 200, a~f = 90, r}„, = 0.3, cr,, 



: 0.05, 



= 45° and a. 
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To investigate the possible range of k and i^o in the two den- 
sity patch model, we have calculated the integrated profiles and 
fitted the FDB curves for 29070 sets of models with different 
combinations of model parameters of a = (5°, 45°, 85°), (3 = 
(1°, 5°, 9°), 7,„ = (200, 300, 400), cr^ = (10, 50, 90), r]„ = (0.1, 
0.3, 0.5, 0.7, 0.9), cr^ = (0.05, 0.1, 0.15, 0.2, 0.25), 4>,„ = (25°, 
45°, 65°, 85°) and ct^ = (5°, 15°, 25°, 35°, 45°). The distri- 
butions of the power-law indices k and the characteristic fre- 
quency i/Q are shown in Fig. [16] We noticed that k values are in 
the range between -0.1 and -2.5. Two components are always 
resolved, except for the extreme case shown in Fig. \T7\ a) in 
which emission is radiated by particles with a narrower energy 
distribution (a^ = 10) and peaking at a larger energy (7„, = 400) 
at very high ri,„ = Q.9 and (j>,„ = 65. 

In above beam calculations, we used the profile width Wio at 
10% of pulse peak. If the separations Wpp of two components 
are used to calculated the beam radii, the beam sizes are ob- 
viously smaller than those calculated from the 50% and 10% 
widths (Fig. [Tsl). We compared in Fig. [19] the k and values 
calculated from i?io (written as kio and I'o.io) with the k and j/q 
values from dpp (written as kpp and v^ pp), and found that they 
are correlated and scatter at small V() or large k values (due to 
error-bar as expected). 

In summary, the curvature radiation from particles in two 
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Fig. 18. — Beam radii (i?io) calculated from the 10% profile widths (Wio) 
at a series of frequencies are compared with beam radii (i?5ci) calculated from 
the 50% profile widths (W50) and beam radii {•dpp) calculated from component 
separations iWpp). Parameters for this example model are 7,,, = 200, = 50, 

= 5°. 



r),„ = 0.7, = 0.2, ( 



= 65° and cta = 15°, with P = Is, a = 45° and /3 





Fig. 19. — Comparison of k and uq values calculated from i?io and 1 



density patches can not only give a proper range of k values, 
but also keeps the two components resolved in general even at 
very high frequencies. 

4. DISCUSSIONS 

In the previous section, we conclude that the FDB curves 
with two resolved components can be preferably explained by 
curvature radiations of particles in two density patches. Here 
we discuss influences on the FDB curves by the radial decay 
of particle energy and by other possible energy distribution of 
particles, and discuss the possible FDB curve flattening at low 
frequencies. Finally we will look at the real FDB data and make 
a model. 

4. 1 . Radial decay of particle energy and the lowest emission 

height 

As shown in Fig. [9] and [15] the values of beam radii at the 
infinite-frequency, i^o, calculated from models are always very 
close to those of /3, no matter what model parameters are in 
putted. However, both on observational ground (Table. [T) and 
theoretical expectations, d{) is larger than (3. 

In simulations above we assume that particles in the mag- 
netosphere have a Gaussian energy distribution with a peak 
Lorentz factor 7,,, and a spread width along the whole tra- 
jectory. In reality, the particle energy distribution^ may be much 
more complicated due to the Compton loss (Zhang et al.ll 19971) 
or other processes. It is reasonable to believe that particle en- 
ergy decreases when they flow out along the field lines in the 
inner magnetosphere. 

In the vacuum-gap model dRuderman & Sutherland! [T975I) . 
electron-positron pairs are produced by the avalanche process 
in the pulsar polar cap. The maximum Lorentz factor of the 
primary particles after accelerating across the gap is 7niax = 
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Fig. 20. — FDB curves for the curvature radiation model b y part icles in two 
density patches with an energy distribution of E g. jg) and Eq. \16\ . The default 
model parameters are the same as those in Fig |13l (e.g. /3 = 3°) but also with 
7„,o = 10*, g = 0.3 and 7., = 350. 



1.2 X 10'^ B^^nhl/Pu, with = B^/IO'^G the surface mag- 
netic field, h4 = /z/lO^^cm the gap height and Pis = P/ls. Pair 
cascade occurs within a few stellar radii near the surface, which 
generates the secondary particles with much lower Lorentz fac- 
tors of a few hundred. Detailed numeric al simulations on th e 
pair cascade process were carried out by iMedin & Lail (1201 Ob . 
in which they found that for various initial parameters (such as 
surface magnetic field, period and primary energy), the pair cas- 
cade process maybe finishes at a height of several star-radius, 
from less than IR^, to more than 10/?*. Here we naively assume 
a toy model to describe the rapid decay of particle energy dur- 
ing the pair cascade process, which reads 



7m = 7mO exp(-^- 



-) + 7.5, 



(16) 



here jm is the peak Lorentz factor of the Gaussian distribution 
we used above, 7^0 and 7^ are the corresponding values of the 
primary particles and the final secondary particles before and 
after the pair cascade, respectively. The factor q (less than 1) 
determines how fast the energy decreases w ith height for the 
cas cade. This en ergy model has been used by Qiao et al] (1200 ll) 
and IZhang et alj (2007) to simulate the inverse-Compton scat- 
tering emission, where they ignored the energy term 75. To use 
the above energy model in our simulation for the model for two 
density patches, one should substitute 7,,, of Eq. ( fT6l ) into the 
Gaussian energy distribution Eq. (|9]l. Note that the number of 
the particles «eo in Eq. ( fTSl l should also be replaced by «eo7.s/7m 
to keep the conservation of total energy. 

Example FDB curves from the curvature radiation model by 
particles in two density patches with an energy distribution de- 
scribed by Eq. (|9]l and Eq. ( fTSI l are demonstrated in Fig. |20] 
and the fitted parameters are listed in Table. 2] We found that a 
larger primary Lorentz factor 7,„o or a smaller q lead to larger 
beam radius at very high frequencies. A larger 7, leads to a 
steeper FDB curve. 



Table 4 

Fitted parameters of FDB curves in Fig.[201 
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Fig. 21. — Characteristic Lorentz factor jc at lowest radio emission height 
well related to ■y„,Q and q. Symbols * come from simulation data as in Table|4] 
and lines are the fitting. 



In these model the values of 1)0 are always larger than those 
of (3, which means that there exists the lowest radio emission 
height rio^v = /?ic sin^(2??o/3), which is > R\cSm^(2/3/3). For a 
small value of -ffo, the corresponding beam radii at the infinite 
frequency can be written as -ffo — 3/2y^now7^- Now, what is 
the value of How? 

According to Eq. ( fTSI l. at the lowest radio emission height 
riow, there exists a characteristic Lorentz factor 7^. We noticed 
that now and jc are related with 7,„o and q by; 

'}^ = l-}.ln^. (17) 

R* q 7mO 

From simulations similarly as presented Table 2] we can get a 
set of riow from do for various 7,„{) and q, and then we got 7c 
using Eq. ( fTTI l. We find that 7c are also very well related to 7„,o 
and q, as shown in Fig. [21] by 

7c = 55.07l4exp(-2^). (18) 
Substituting into Eq. ( fTTT i. we find 



L381og7^ 
How =/?*(3 + ). 

q 



(19) 



4.2. On the energy distribution of particles 

At any height, the secondary particles in the magnetosphere 
for the curvature radio emission were assumed to have an en- 
ergy distribution in the Gaussian function. It is possible that the 
particles follow other energy d istribution functions. For exam- 
ple, according to lAronsl(ll98lb . the multi-component plasma in 
pulsar magnetosphere contains secondary pairs, a high-energy 
plasma 'tail', and some primary particles. The joint energy 
distribution of secondary particles and high-energy 'tail' can 
be described by a power-law with two cutoffs at its two ends, 
which reads as, 

«e(7) ~ 7''7 7i,min < 7 < 7i,max, (20) 

where s is the power-law index. The curvature radiation from 
particles in two density patches with such an energy distribu- 
tion produces the integrated profiles at a series of frequencies 
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Fig. 22. — Same as Fig. ll4l but a power-law energy distribution with index 
i = -2. Other model parameters are 7, niin = 150 and 7i,max = 2000, 7,,, = 350, 
CT^ = 30, »?,„ = 0.75, o-^ = 0.08, (/i,„ = 45° ando-^ = 12.5°, withP= Is, a = 45° 
and ^ = 3°. 
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Fig. 23. — Example of Low frequency flattening the FDB curve. The model 
parameters are 7™ = 400, = 10, rj,,, = 0.7, <Jrf = 0.1, (/>„, = 45° and cr^ = 5° 
with P = li, a = 45° and /3 = 5°. 



with a similar shape and a constant pulse width between 80 
MHz and 40 GHz, as shown in Fig. |22l We can not find the 
profile differences for a range of s values, e.g. a difference less 
than 2% for the pulse width values between s = -2 and s = -3. 
We also found that the values of 7j.inin and 'js.max determine the 
frequency range for the constant profile width. 

4.3. Possible low frequency flattening of the FDB curves 

During the model calculations, we noticed that the FDB curves 
tend to be flattening at low frequencies. For some input param- 
eters of the two density patch model, the effect is very obvious 
(Fig.l23]l. especially when Lorentz factors are large (7,,, > 400) 
and/or the energy distribution is broad {a^ > 80). 

This is understandable. Because of the dipole bending of 
magnetic field lines, as pulsar rotates, a fixed sight-line cut 
across the two density patches (see Fig. [TZt with a clear lower 
limit and upper limit of height and a clear left and right limit 
of the rotation phase. Very low frequency emissions tend to 
come from the high regions of two patches in such a limited 
volume. In that high regions near the upper bound, the lower 
frequency emission comes from particles of the lower energy, 
but the pulse-width or beam-width of the emission from two 
patches is determined by the geometry. This happens to very 
low-frequency, depending on the 7,„. The low-frequency flat- 
tening of FDB curves are often seen in the models when patches 
are small (small cr^), or energy of particles are very high and/or 
the energy distribution is broad (large 7^ and a~^). 

4.4. Observed FDB curve and a model 
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Fig. 24. — Upper panel: An example of rea l data and model - fitting . The 
observed data of PSR B2020-I-28 are taken from lMitra & Rankinl jlOOa . The 
solid line is the best model calculated from the curvature radiation of particles 
in two density patches with parameters 7,„o = 10^, 7.1 = 225, = \5,q = 0.25, 
r},n = 0.75, (T^ = 0.08, <!>„, = 45° and = 12.5°. The geometrical parameters 
a and /3 and the period of this pulsar can be found from Table 1 . Lower panel: 
The X' distribution for a pair of the most influential parameters, 7, and q, for 
the fitting. The best values are 7j = 225 and q = 0.25. 



We have worked on a huge number of models for the radio 
emission intensity distribution across the height and rotation 
phase and integrated profiles. We now look at data of a real 
cone-dominant pulsa r, PSR B2020-H28, and make a model. The 
data are taken from (iMitra & Rankinl |2002|) . To keep the two 
profile components separated even at very high frequencies, we 
have to choose a model of the curvature radiations of particles 
in two density patches and consider the energy distribution of 
particles as described by Eq. (|9) and Eq. (fT&t . As shown in 
Figure |24] the model can fit data nicely with a local minimum 
of = 10.24 in the parameter space of 7, and q. The model 
is obtained with sparsely separated parameter grids, which may 
be improved by a global fitting for 8 parameters but that is very 
computational expensive. 

5. SUMMARY AND OUTLOOK 

We studied the frequency dependence of beam radius for 
cone-dominant pulsars by using the curvature radiation mecha- 
nism. 

For the simplest case in which pulsar radio emission is gener- 
ated by curvature radiation of relativistic particles with a single 
7 and streaming along fixed open field-lines, we obtain the an- 
alytic formula in Eq. d?) for the frequency dependence. It has 
the power-law index k = -\, and i?o = 0, both are different from 
observations. 
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Considering various density distribution and energy distri- 
bution of particles, we numerically calculated the intensity of 
curvature radiation from all the possible emission heights in the 
open field line region, and get integrated mean pulse profiles 
from which the FDB curves are derived. We found that the 
power-law index k and are strongly influenced by the den- 
sity distribution and energy distribution of particles in the open 
field-line region. 

When particles are uniformly distributed in whole open field 
lines, the curvature radiation can produce only one profile com- 
ponent, which can not explain the two resolved components of 
cone-dominant pulsars. 

When particles are distributed in a conal area, the low-frequency 
emission comes from very high regions so that profile is very 
wide. Profiles at high frequencies cannot have two separated 
components, because emission comes from low heights where 
the sight line always continually cuts the edge of density cone. 

To keep two profile components separated or resolved, the 
particles must come from two density patches. Model calcu- 
lations from a huge amount of parameter combinations show 
that the k values of the FDB curves are in the range of -0. 1 to 
-2.5, in agreement with the observations. Because of the lim- 
ited upper boundary of the density patches cut by a given sight 
line, the FDB curves tends to be flattened at low frequencies, 
especially for models with a large Lorentz factor and/or broad 
particle energy distribution. 

The geometry of density distribution of particles and the en- 
ergy distribution of particles in the model certainly affect the 
FDB curves. For a Gaussian energy distribution, high energy 
particles with large 7^ will produce a wide beam at low fre- 
quency. Particles in a narrow energy distribution with a small 
cr-y produce a steep FDB curve. The radial decay of particle 
energy is important to get a reasonable ?9o in models. 

Through model calculations of pulsar beam and profile evo- 
lution by using the curvature radiation mechanism, we have 
clarified some important issues on particle density and energy 
distributions. However, in our simulations, we considered only 
the simple case that particles are emitting radio waves with the 
characteristic frequency for curvature radiation determined by 
the Lorentz factor. The real curvature emission of any particles 
has an energy spectrum, which is a power-law with an index 
1/3 at the low frequency limit (much less than the character- 
istic frequency) and becomes exponential at high frequencies. 
This may significantiy affect the emission regions and modify 
the FDB curves, which we should investigate in future. Another 
important fact is that the real pulsar radio emission is believed 
to be highly coherent, while in this paper we simply treat the 
coherence by assuming that the emission comes from point-like 
huge charge. Moreover, other emission mechanisms may also 
work in pulsar magnetosphere, such as inverse Compton scat- 
tering, plasma oscillation, cyclotron instability and Coulomb 
bremsstrahlung emission, etc, which are not considered in this 
paper and may lead to their own FDB curves. 

The authors thank Professor Chou Chih Kang for helpful dis- 
cussions and the referee for constructive comments. This work 
has been supported by the National Natural Science Foundation 
of China (11003023 and 10833003). 
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